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Abstract. We present results of modelling the bulk of the spectral energy distribution (0.35 - 5 fim ) for CJ406 
(M6V). Synthetic spectra were calculated using the NextCen Dusty and Cond model atmospheres and incorporate 
line lists for H2O, TiO, CrH, FeH, CO, MgH molecules as well as the VALD line list of atomic lines. A comparison 
of synthetic and observed spectra gives T^s — 2800 ± 100 K. We determine Mboi = 12.13 ± 0.10 for which 
evolutionary models by Baraffe et al. (2003) suggest an age of around 0.1 - 0.35 Gyr consistent with its high 
activity. The age and luminosity of GJ406 correspond to a wide range of plausible masses (0.07 - 0.1 Mq ). 
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1. Introduction 

Studies of M dwarf spectra are of interest to many 
. branches of modern astrophysics. Indeed, perhaps 70 % 
of stars within 10 parsecs are M dwarfs and it is very 
. probable that this number density prevails throughout our 
I Galaxy. The population of these numerous low-mass stars 
! (0.08 Mq < M < 0.6 Mq), together with substellar ob- 
] jects (brown dwarfs; M < 0.075 Mq ) would contain an 
. appreciable amount of the baryonic matter in the Galaxy. 
' Estimates of brown dwarf number densities currently sug- 
. gest the same order as for stars {^^ 0.1 per pc~'^), therefore 
' their contribution to the total mass should not exceed 15 
1 % (Reid et al. 1999). Nonetheless, the large errors associ- 
ated with age and mass determinations for brown dwarfs 
make such estimates very uncertain. 

The verification of the theory of stellar evolution and 
structure of stars, the detection among M dwarfs of a 
subset of young brown dwarfs, and the physical state of 
plasma in their low temperature atmospheres are among 
a few of the interesting problems that may be addressed 
through the detailed study of M-dwarfs. 

Some authors reference GJ406 (other names are V*CN 
Leo, EUVE J1056-H07.0, [GKL99] 228, GSC 00261-00377, 
LET 750, LHS 36, 2MASS J10562886-H0700527, 2RE 
J1056-I-070, IRXS J105630.3-h070118) as an "archetype 
dwarf of spectral type M6V" , or one of the "well known 
spectral standards for its type" (Mohanty et al. 2004). 

Send offprint requests to: Ya. V. Pavlenko 
Correspondence to: yp@star.herts.ac.uk 



GJ 406 is located at 2.39 pc from the Sun (Henry et al. 
2004). Altena et al. (1995) determined a proper motion of 
fj, = 4.696. Radial velocities are of order 19 ± 0.1 km/s 
(Martin et al. 1997, Mohanty & Basri 2003, Euhrmeister 
et al. 2005). Leggett (1992) found that this nearby dwarf 
has typical old disk properties. Deflosse et al. (1998) re- 
ported a rather low v * sinz < 3 km/s (see also Mohanty 
& Basri 2003). Guetter et al. (2003) even use GJ 406 as 
one of the JHK standard stars on the CIT photometrical 
system. Indeed, optical and IR spectra of the dwarf do not 
contain any unusual features. They are governed by ab- 
sorption of diatomic molecular band systems, such as TiO 
and VO (Kirkpatrick, Henry & McCarthy 1991), as weU 
as rotational- vibrational bands of H2O and CO (Jones et 
al. 1994). 

On the other hand, GJ 406 has 
some very unusual properties. SIMBAD 
(http://simbad.u-strasbg.fr/sim-fid.pl) labels GJ 406 
as a flare star. The dwarf has strong H^ of EW — 6.7 
A(LH^/Lboi = -3.9, see Mohanty & Basri 2003). An 
ultraviolet spectrum of GJ406 contains some emission 
lines (Euhrmeister et al. 2004). Furthermore, Schmitt & 
Wichmann (2001) detected the Fe XIII forbidden coronal 
line at 0.33881 fim . Euhrmeister at al. (2004) reported 
a high level of variability of this line on a timescale of 
hours which they ascribe to microflare heating. Recently 
an X-ray luminosity of log = 26.97 was detected 
by Schmitt & Liefke (2004). GJ 406 is the only known 
M6 star yet observed with a strong chromospheric and 
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coronal activity. Only a few stars in the solar vicinity are 
known with such a menagerie of activity phenomena. 

In this paper we compute the synthetic energy distri- 
bution of several model atmospheres with a range of effec- 
tive temperature and compared them with the observed 
fluxes of GJ406. Section 2 presents the spectral data used 
in our paper. Section 3 describes our procedure for com- 
putation. We psesent our results in Section 4. In Section 
5 we discuss the implications of our results. 

2. Observations 

Table m lists the data and instruments used to obtain the 
observational spectra in this paper. The spectra are shown 
in Fig. n and come from measurements taken with a va- 
riety of different instruments on different telescopes. All 
the data are assessed to be of good quality, most have 
already be used for other papers. We refer to the spec- 
tra in wavelength order. The reduction procedures for the 
Integral Spectrometer (ISIS) data taken on the William 
Herschel Telescope (WHT) are described in Dobbie et al. 
(2004). The reduction procedures for the CGS4 data taken 
on the United Kingdom Infrared Telescope (UKIRT) here 
have been reported elsewhere in Jones et al. (1994, 1995, 
1996). It should be noted that CGS4 observations in the 
spectral regions labelled in Fig. by "m" and "n" are not 
continuous. Fluxes between them were filled by NextGen 
synthetic spectra. The reduction procedures for the Short 
Wavelength spectrometer (SWS) data taken with the 
Infrared Space Observatory (ISO) are reported elsewhere 
in Jones et al. (2002). The Near Infrared Camera Mosiac 
Spectrograph (NICMOS) data taken on the Hubble Space 
Telescope (HST) was reduced using the data processing 
software Calnicc (Version 2.5.7). The data has been com- 
pared with Jones et al. (1994) and is preferred due to 
its excellent flux calibration (5-10%, Pirzkal & Freudling 
1998) and broad wavelength coverage. 

In general we have renormalised the fluxes for the dif- 
ferent spectral regions so as to ensure that flux levels are 
the same where the regions overlap. For the CGS4 spec- 
tra centered on 4.7 /xm (region n), this was not possible 
since there is no overlap. For this spectral region, we sim- 
ply used the instrument calibration to determine the flux 
level, and filled in the gap with a normalised section of 
NextGen model spectra. 

2.1. Absolute flux calibration. 

In order to provide an absolute normalisation for the full 
GJ406 spectrum, we used the available near infrared pho- 
tometry from Leggett et al. (2000), which we transformed 
onto the Mauna Kea Observatory (MKO) photometric 
system using Hawarden et al. (2001). We used the mea- 
sured MKO J, H and K filter pass bands convolved with 
atmospheric transmission to estimate the ground based 
photometrically measured flux components in our spec- 
tra, by integrating over each band. We then performed the 
same task on a flux calibrated spectrum of Vega. Vega was 



assumed to be zero magnitude at all wavelengths, and the 
GJ406 flux could thus be scaled to match the photometry. 
The normalisation value that we found was ^20% higher 
when normalising in the K-band than in the J-band, and 
was an intermediate value in the H-band. These differ- 
ences presumably result primarily from the relative nor- 
malisations that we used to join the individual spectral 
regions, and provide an accuracy gauge on this procedure. 
We chose to make our final normalisation in the H-band, 
and duly estimate a likely ±10% uncertainty in our abso- 
lute flux levels. 

In order to derive a bolometric flux, we added a 4.8- 
20 /im spectral tail to our calibrated spectrum (using 
a NextGen 2800 K, logg=5.0, [M/H]=0 model spectra, 
shown to be appropriate in Section 4.2), normalised in its 
overlap region, and integrated out to 20 /im . Note how- 
ever, that this synthetic spectral tail only contributes ^3% 
to the bolometric flux, which we found to be 6.35 x 10^^^ 
W tot'^. We then derived the bolometric magnitude (mboi) 
using the Sun as a standard (adopting 7^0=3.86 x 10^^ W 
and Mboi0=4.75), which yields niboi=9.02. This is consis- 
tent with the value derived by Leggett et al. (2000) of 9.07. 
Assuming a distance modulus of m = -3.11 ± 0.01 (van 
Altena et al. 1995) and 10% uncertainty in our flux calibra- 
tion, we thus determine that GJ406 has Mf,oi=:12.13±0.10 
and L = logL*/Lo=-2.95±0.05. 

3. Theoretical spectra computation procedure. 

Theoretical spectral energy distributions^ were computed 
for model atmospheres of dwarfs with effective temper- 
atures Toff = 2500-3200 K from the NextGen grid of 
Hauschildt, AUard & Baron (1999) for solar metallic- 
ity (Anders & Grevesse 1989). Hereafter we use the 
syntax "effective temperature/ gravity /metallicity", e.g 
2800/5.0/0 to signify the model atmosphere. Unless 
otherwise mentioned all models are for log g — 5.0. 
Computations of synthetic spectra were carried out by 
the program WITA6 (Pavlenko 2000) assuming LTE, hy- 
drostatic equilibrium for a one-dimensional model atmo- 
sphere and without sources and sinks of energy. The equa- 
tions of ionisation-dissociation equilibrium were solved for 
media consisting of atoms, ions and molecules. We took 
into account ~ 100 components (Pavlenko 2000). The con- 
stants for equations of chemical balance were taken from 
Tsuji (1973). 

Molecular line data were taken from different sources. 
Lines of ^Hj^O were computed using the AMES database 
(Partrige & Schwenke 1998). The partition functions of 
H2O were also computed from these data (see section 
3.1). i^ci^o and ^^C^^O hue lists were computed by 
Goorvitch (1994). The CO partition functions were taken 
from Gurvitz et al. (1989). TiO line lists were taken from 
Plez (1998) and Schwenke (1998). CN lines came from 
CDROM 18 (Kurucz 1993); CrH and FeH lines were taken 

^ Hereafter we use the term "synthetic spectra" to simplify 
the text. 



Pavlenko et al.: SEDs of GJ406 3 
Table 1. Telescope and Instrument configurations used to collect our GJ 406 dataset. 



Start End Instrument (configuration) Telescope Date 
jum jj.m 



0.35 





56 


ISIS (blue arm) 


WHT 


2001 


Jan 29 


0.55 





80 


ISIS (red arm) 


WHT 


2001 


Jan 29 


0.79 


1 


20 


NICMOS fG096) 


HST 


1998 


June 19 


1.05 


1 


95 


NICMOS (G141) 


HST 


1998 


June 19 


1.3 


2 


59 


NICMOS (G206) 


HST 


1998 


June 19 


2.48 


2 


60 


SWS (06 1A~) 


ISO 


1996 


June 26 


2.60 


2 


75 


SWS (06 lA) 


ISO 


1996 


June 26 


2.74 


2 


90 


SWS (06 lA) 


ISO 


1996 


June 26 


2.88 


3 


02 


SWS (06 IB) 


ISO 


1996 


June 26 


3.03 


3 


23 


CGS4 (150 1/mm) 


UKIRT 


1993 


April 20 


3.21 


3 


40 


CGS4 (150 1/mm) 


UKIRT 


1993 


April 20 


3.36 


3 


75 


CGS4 (75 1/mm) 


UKIRT 


1992 


May 7 


3.76 


4 


15 


CGS4 (75 1/mm) 


UKIRT 


1992 


May 7 


4.51 


4 


90 


CGS4 (75 1/mm) 


UKIRT 


1992 


October 26 



7 
6 
5 



X 

LJ- 3 



2 - 



T r 




a) ISIS (blue arm) on WHT at 2001 Jan 29 
b) ISIS (red arm) on WHT at 2001 Jan 29 

c) NICMOS (G096) on HST at 1 998 June 1 9 

d) NICMOS (G1 41 ) on HST at 1 998 June 1 9 

e) NICMOS (G206) on HST at 1 998 June 1 9 

f) SWS (06 1A) on ISO at 1996 June 26 

g) ibid 

h) ibid 

i) SWS (06 1 B) on ISO at 1 996 June 26 
j) CGS4 (150 1/mm) on UKIRT at 1993 20 April 

k) ibid 

I) CGS4 (75 1/mm) at 1992 May 7 

m) ibid 

n) CGS4 (75 1/mm) 1992 October 26 




0.5 1 1.5 2 2.5 3 3.5 

Wavelength (micron) 



4.5 



Fig. 1. Observational data used for this paper. The wavelength coverage of the different instruments is shown. 



from Burrows et al. (2002) and Dulick et al.(2003), respec- 
tively. Atomic line list was taken from VALD (Kupka et 
al. 1999). 

The profiles of molecular and atomic lines were de- 
termined using the Voigt function H{a,v). Parameters of 
their natural broadening C2 and van der Waals broaden- 
ing C4 were taken from Kupka et al. (1999) or in their ab- 
sence computed following Unsold (1955). Owing to the low 
temperatures in M dwarf atmospheres and consequently, 
electron densities, Stark broadening could be neglected. 
As a whole the effects of pressure broadening prevail. 
Computations for synthetic spectra were carried out with 



a step 0.5 A for microturbulent velocity ut = 1 - 4 km/s. 
The instrumental broadening was modelled by gaussian 
profiles set to approximate the resolution of the observed 
spectra. The relative importance of the different opacities 
contributing to our synthetic spectra is shown in Fig. El 

3.1. Partition functions of water 

We recomputed the constants of chemical equilibrium fol- 
lowing Kurucz (1970) taking into account weights Si = 
1/4 and 3/4 for levels of water of different symmetry. We 
followed the scheme described by Pavlenko (2002). Let us 
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Fig. 2. The contribution of different molecules to the 
formation of the synthetic spectrum of in Tcff /log g = 
2800/5.0 model atmosphere. 



write an equation of ionisation-dissotiation equilibrium for 
the molecule consisting oi x, y,...., z atoms as 



^n^jn^,,,^ = exp{-E.^y,,,^/Tev + h- c* 



iT + d*{T^e^{T + f*T))) + 
3/2 * (to - A: - 1) * InT) 



(1) 



where E and T are dissociation energy and temperature 
(in eV), Uz is the number density of z-species, k and m 
are ionisation degree (0 for neutrals) and number of atoms 
per molecule, respectively (see Kurucz (1970) for more de- 
tails). Computed constants a, 6, c, d, e, / are given in Table 

m 

The temperature dependence of computed partition 
functions of H2O are given in Fig. |3 In general, our 
new partition functions agree well with data of Vidler & 
Tennyson (2000) computed using a mixture of experimen- 
tal data and a UCL model of the water vapour molecule. 
Some differences occur at T > 5000 K as the UCL model 
has more levels of high excitation energy. 
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Fig. 3. Partition functions of molecules H2O and 
HDO. Our data for H2O are compared with Vidler & 
Tennyson(2000), data for HDO are compared with Hewitt 
et al. (2005). Our partition function for HDO is computed 
for Si =1. The differences in the HDO partition functions 
at higher temperatures (T > 4000 K) are due to the use of 
more complete sets of deuterated water levels in the UCL 
model compared with the AMES model. 
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Fig. 4. Molecular densities of H2O computed for NextGen 
model atmospheres 2800/5.0 and 3800/5.0 for the cases of 
partition functions computed with Si = 1,4 and Si =1/4, 
3/4. 



It is worth noting that 
- the molecular densities of given molecules obtained from 
a solution of the system of equations of molecular equilib- 
rium response to changes of absolute values of U, espe- 
cially in high temperature regimes (T > 3000 K, see Fig. 
HI). 

- J7(H2 ) depends strongly on temperature. Table [21 
provides fitting constants obtained for different temper- 
ature regions. From a general point of view it would be 
reasonable to restrict our fitting to temperatures T > 300 
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Table 2. Constants of the dissociation equilibrium of H2O in formats of the ATLASxx, where xx labels a version of 
ASTLAS, the superscripts correspond to different fit temperature ranges: ^ — 300 K < T < 1000 K, ^ — 300 K < T 
< 6000 K, 3 — 60 K < T < 6000 K, * — 50 K < T < 10000. 



Do 


h 


c 


d 


e 


/ 


Refs. 


9.500 


0.9193E+02 


0.2550E-02 


0.4088E-06 


0.3893E-10 


0.1512E-14 


Pavlenko (2002)^ 


9.512 


9.3179E+01 


2.6725E-03 


5.7830E-07 


8.5268E-11 


5.1311E-15 


Kurucz (1999) 


9.500 


0.9374E+02 


0.3494E-02 


0.9795E-06 


0.1655E-09 


0.1073E-13 


Vidler & Tennyson (2000)^ 


9.500 


0.9494E+02 


0.5858E-02 


0.2338E-05 


0.4634E-09 


0.3279E-13 


Vidler & Tennyson (2000f 


9.500 


0.9331E+02 


0.2551E-02 


0.4089E-06 


0.3895E-10 


0.1513E-14 


This work^ 


9.500 


0.9428E+02 


0.3722E-02 


0.8305E-06 


0.9719E-10 


0.4238E-14 


This work 



K. Though for some astrophysical objects it would be in- 
teresting to have the partition function of water vapour 
for even lower temperatures. Thus we have provided these 
as well. 

4. Results 

4.1. Dependence of theoretical spectra on input 
parameters 

First of all, we computed model SED's using different in- 
put parameters: effective temperatures, gravities, metallic- 
ities, microturbulent velocities. Ratios of fluxes computed 
with different sets of input parameters are shown in Fig. 
O It can be seen that a temperature change of 200 K is 
roughly equivalent to a change in metallicity of 0.5 dex or 
a gravity change of Alog g — 1. 

In Fig. we see the differential effects of the depen- 
dence of our model spectra on different parameters. Some 
of these effects depend more on changes in the opacities for 
different parametres. However, some effects are explained 
by changes in the structures of the model atmosphere. 
Indeed, stellar photospheres of different Tes , log g, [M/H] 
lie in different pressure regions. 

4.2. Dependence of theoretical spectra on different 
model atmospheres 

We also examined the dependence of our results on the 
choice of model atmospheres: NextGen, Dusty, Cond (see 
AUard 2005 for references) . These model atmospheres have 
different T = f{P) structures (Fig.|HJ due to differences in 
the physical treatment of the dust formation which cause 
some changes in the opacities and the molecular equilib- 
rium. In general, the DUSTY and COND models have 
more hot inner layers and cooler outer layers (see first 
plot in Fig. EJ. 

We computed model theoretical fluxes for the 
NextGen, Dusty and COND model atmospheres with pa- 
rameters 2800/5.0/0 in the spectral region of interest and 
compared them. Theoretical fluxes were convolved with a 
spectral resolution element of FIIWM=4 A. In Fig. we 
show the ratio of the convolved fluxes. 
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Fig. 6. Top: Temperature structure of model NextGen, 
DUSTY and COND model atmospheres of 2800/5.0/0. 
Bottom: ratio of fluxes computed for these model atmo- 
spheres. 

It is worth noting a few results: 

— Infrared spectra containing strong water bands 
agree rather well. 

— We see rather big differences in the optical and blue 
parts of the spectra. Here the dependence of computed 
spectra to changing temperature is much stronger in com- 
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Fig. 5. Responses of computed spectra to variations of input parameters. A model atmosphere of 3000/5.0/0.0 was 
used as the reference model. 
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Fig. 7. Ratio of fluxes computed for 2800/5.0 NextGen 
model atmosphere with TiO line lists by Plez (1998) and 
Schwenke(1998). 



parison to the IR spectral region. In general, both the 
COND and DUSTY model atmospheres are a bit hotter in 
the inner regions, and we therefore see an increased flux in 
the continuum flux levels (or increased flux for lower opac- 
ity wavelength ranges). Some molecular lines here become 
stronger due to the lower temperatures in DUSTY and 
COND model atmospheres. 

— CO bands strength responds to the change in tem- 
perature in the photospheric layers, due to the high sensi- 



Fig. 8. Min F(/s,/„) found for different model atmo- 
spheres. 

tivity of dissociation equilibrium of this molecule to tem- 
perature. 

4.3. Dependence of theoretical spectra on the use if 
different TiO line lists. 

There arc two TiO line lists of widespread use by 
Schwenke(1998) and Plez(1998). They are based on the 
improved Langhofl[(1997) model of the TiO molecule but 
difii'er in details. Plez (1998) added an a-f system at 0.5 
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Fig. 9. Fits with linear and logarithmic flux scales to our GJ406 spectrum with a theoretical spectrum computed for 
a solar composition NextGen model atmosphere 2800/5.0. 



/zm to the line list. Schwenke subsequently computed a 
corresponding list of transitions complete to the higher 
excitation energies. Plez (1998) provided a line list for a 
solar mixture of Ti isotopes. Schwenke's (1998) database 
provides lists for each TiO isotope. 



We compare synthetic spectra computed with these 
two line lists for a NextGen model atmosphere 2800/5.0/0 
(Fig. d). Synthetic spectra were computed with a 0.1 
A step, then convolved with a spectral resolution element 
of FWHM = 1 A. 



We find the largest differences in the blue part of 
spectrum, which is more affected by incompleteness of 
molecular line lists for other molecules, chromospheric 
effects, veiling, and strong atomic absorption. Therefore 
from inspection of Fig. [7|we conclude that differences be- 
tween synthetic spectra computed with Plez (1998) and 
Schwenke (1998) line lists do not affect our main results 
(see also Lyubchik & Pavlenko 2001). 



4.4. Fits to GJ406 spectra 

In this paper we are interested in the dependence of com- 
puted spectral energy distributions to T^g . To find the 
best fit of computed spectra to observed fluxes we use a 
minimisation procedure described in Jones et al (2002) 
and Pavlenko & Jones (2003). Namely, the best fits are 
found for the min F{fs,fn,fg), where fs, fn and fg are 
relative shift of the spectra, normalisation constant for 
the computed spectra, and broadening parameter, respec- 
tively. We found a rather weak dependence of the min F 
on the fg broadening parameter and set /g = 6 A. Previous 
studies have considered GJ406 as a typical M6 dwarf, and 
we thus assumed log g = 5.0 in its atmosphere. 

To determine a self-consistent solution we fit the theo- 
retical spectra to the observed fluxes in all spectral regions 
and estimate the quality of the flt by computing F{fs, /„). 
Two spectral regions were excluded from our analysis: 0.35 
- 0.4 /im due to incompleteness of our opacity sources, and 
4.3 - 0.461 /im due to the gap in the observed data. 

In Fig. IHl we show computed F(/s,/„) for a grid of 
our theoretical spectra of different T^s ■ We find a weak 
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Fig. 10. A comparison of f{xs,Xn) computed for the fit- 
ting of onserved spectra by NextGen and Cond synthetic 
spectra. 



dependence of F{fs, fn) on Vt. The best fit can be found 
for the min /„) at 2800 ± 100 K ( Fig.lHlshows these 
fits in Unear and logarithmic scales). 

Some problems with fitting spectral features are seen 
at 1.3 < A < 1.7 /zm . Partially these discrepancies can be 
exlained by problems with modelling of strong H2O bands 
located here (Jones at al. 2005). Fits to TiO and VO in 
optical spectral regions and H2O (beyond 1.7 /xm ) are of 
particularly good quality (see Pavlenko 1998, Jones et al. 
2002 for more details). 

We also made comparison of the observed spectrum 
with synthetic spectra based on the Cond model atmo- 
spheres and obtained the fits practically of the the same 
quality. The "best fit" synthetical spectra computed with 
NextGen and Cond model atmospheres coincide over the 
whole spectral region. The largest differences do not ex- 
ceed 1-2 %. Consequently, the function f{xs,Xn) for fits 
to the Cond synthetic spectra has the minimum at 2800 
K, as for NextGen (see Fig. lTUIl . 

To confirm this explanation of the discrepancy between 
theory and observation, it would be desirable to carry out 
a similar analysis of late- type dwarfs like GJ 406, but em- 
ploying a more complete water vapor line list for the theo- 
retical atmosphere calculations. Such a line list is expected 
to be available in the near future (Jones et al. 2005). 

4.5. Evolutionary model fits 

Using evolutionary models by Baraffe et al (2003), we have 
estimated a mass and age for GJ 406 (Fig. ITT)). Our mea- 
sured effective temperature Tgff =2800±100 K, and lumi- 
nosity L = log L^/Lq = -2.95±0.05 agree well with esti- 
mates of other authors (Leggett et al. 2002; Jones et al. 
2002). We know from the lack of lithium in GJ 406 that its 
age must be > lO^yrs (Magazzu, Martin & Rebolo 1993). 
Indeed, within the error bars we find that the measured 
Tcff and luminosity are consistent with the Baraffe et al 



(2003) evolutionary models for stars of T^ff = 2700 - 2900 
K and L = -3.0, which complies an age of 0.1 - 0.3 Gyr, 
and, respectively, M = 0.07 - 0.94 Mq. For L = -2.9 we 
obtain the age 0.2 - 0.35 Gyr, and M = 0.1 - 0.09 Mq. 

5. Discussion 

A comparison of observed and computed spectral energy 
distributions provides a unique tool with which to assess 
the completeness and quality of our knowledge about the 
structure and properties of late type dwarf atmospheres; 
the physical state of their matter; opacity sources and line 
lists, atmospheric temperature and pressure structure, ef- 
fective temperature scales. 

Here we have modelled the spectral energy distribu- 
tion of the M6 dwarf GJ 406 from 0.4-4.9 /xm . The op- 
tical spectrum is formed primarily by absorption of the 
saturated bands of VO and TiO. In general, the response 
of optical fluxes to the variations of input physical param- 
eters is stronger than at infrared wavelengths (see FigEJ. 
However, the infrared spectrum is more sensitive to the 
H2O absorption bands, with fluxes coming from deeper 
atmospheric layers. Therefore, if one is to understand both 
the outer and inner atmospheric structure of cool dwarfs, 
one must simultaneously account for both the optical and 
infrared spectra. In general, we achieved good agreement 
between our theoretical spectra and observation, although 
we note some problems with fits at certain wavelengths 
(0.35 - 0.4 iim. , 1.3 - 1.7 /xm , these are presumably a 
consequence of missing molecular opacities). Our fits sup- 
port the idea that the TiO/VO and H2O line lists covering 
these wavelength ranges are of good quality. Most prob- 
ably, our problems in the blue part of the spectrum can 
be solved with proper fits to strong atomic lines located 
there. Then, the blue part of the spectrum should be more 
affected by chromospheric like phenomena. The detailed 
analysis of these and related problems is beyond the scope 
of this paper. We plan to consider those in forthcoming 
papers. 

Our best fit age range is consistent with the age con- 
straints from depleted lithium. However, GJ406 has been 
kinematically classified as an old disk star, and as such, 
one might expect its age to be greater than ^^600 Myrs 
(the age of the Hyades, which traditionally represents an 
upper age limit for the young disk population; eg. Leggett 
1992). Our model fit age range thus suggest that GJ406 is 
more youthful than this, and that its old disk classification 
should be interpreted solely as a kinematic description, 
and not as an age constraint. This is not contradictory of 
course, since the dispersion in the kinematics of the young 
disk population naturally places some young stars outside 
of the canonical young disk UVW kinematic region. 

Although the spectrum of GJ406 shows strong H^ 
emission and activity (as summarised in Section 1), it 
should be noted that one cannot use this to place any 
strong constraints on age. As Gizis et al (2002) ex- 
plain, one expects a significant spread in the H^ emission 
strength of a population if the age is less than some value 



Pavlenko et al: SEDs of GJ406 



9 




Fig. 11. Top: location of GJ 406 in respect to the Baraffe et al. (1998) evolutionary tracks and Toff = 2800, 2600 
K and L = L^ILq =-2.90,-3.00 planes. The planes shown on plot determine the cuboid of our errors of Teg and L. 
Bottom: location of GJ 406 in respect to evolutionary tracks in coordinates -L, t, M. Tthe thick lines shown on plot 
determine the cuboid of our errors of age, luminosity L^/Lq and mass M*. 



that depends on the stellar colour or spectral type (see 
their figure 11). For GJ406 {V-I = 4.06), this age upper 
limit is greater than the age of the disk, and a high level 
of Ha emission may thus be expected for some M6 dwarfs 
spanning the full age of the disk. It is at least clear that 
the emission properties of GJ406 are not inconsistent with 
our age estimate. 
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